Magnetic and orbital blocking in Ni nanocontacts 
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We address the fundamental question of whether magneto-resistance (MR) of atomic-sized con- 
tacts of Nickel is very large because of the formation of a domain wall (DW) at the neck. Using 
ab initio transport calculations we find that, as in the case of non-magnetic electrodes, transport 
in Ni nanocontacts depends very much on the orbital nature of the electrons. Our results are in 
agreement with several experiments in the average value of the conductance. On the other hand, 
contrary to existing claims, DW scattering does not account for large MR in Ni nanocontacts. 

PACS numbers: 



The strong sensitivity of the current flow between 
two ferromagnetic metals (FM's) separated by a non- 
magnetic region to the relative orientation of their 
magnetization vectors is a fundammental physical phe- 
nomenon with a huge impact in the magneto-electronics 
industry- . The figure of merit is the ratio between the 
conductance for antiparallel (AP) and parallel (P) rela- 
tive magnetic orientations of the FM's, x = ^^, which 
can be selected with an external magnetic field. Two dif- 
ferent conventions are used to characterize the so called 
magneto-resistance (MR), MRi = 100 x (1 — x)% and 
MR2 = 100 X (i - 1)%. Metallic multilayers with al- 
ternated ferromagnetic and non-magnetic metallic layers 
display a large MR known as Giant MR (GMR)^. MR 
in FM/insulator/FM systems is known as Tunnel MR 
(TMR)1 and values of MRi ~ 30% have been reported!. 

More recently a number of gi-oupsS*SiL2i2iiSiii have 
studied MR in break junction systems, where two sec- 
tions of a Ni wire are connected through an atomic-size 
contact. In this arrangement the intermediate region 
connecting the two bulk FM's has a different geome- 
try but the same chemical composition, in contrast to 
GMR and TMR systems. Some groups have obtained 
values of MR2 going from 10^ to 10'* (x ~ 10~^ to 
X ~ j^g-3^^5,6^ while others obtain moderate or even 
negative valuea^i^iiSiii. In nanocontacts the resistance 
predominantly comes from the region with the smallest 
section, where electron transport is coherent and conduc- 
tance is dominated by the quantum mechanical transmis- 
sion of electrons at the Fermi leveUl. Mainly, two differ- 
ent mechanisms have been proposed so far to account 
for the large values of MR, when observed: Domain-wall 
(DW) scatteringi^ii^ and magnetostriction*"'^^. For the 
former it has been argued that in the AP arrangement 
a DW is pinned at the atomic-size contact*^ and is re- 
sponsible for strong spin scattering which gives an extra 
contribution to the resistance as compared to the P con- 
figuration, resulting in a large "ballistic" MRliW-. 

The fundamental question of whether MR is dramati- 
cally enhanced in atomic sized ferromanetic contacts due 
to the presence of a DW remains open and is the sub- 
ject of this paper. Three ingredients are essential to an- 
swer this question. First, as in the case of non-magnetic 
nano-contacts, the electronic structure of the last atom(s) 



which determines the number of available transport 
channels^. Second, the presence of inhomogeneous mag- 
netization profiles, e.g., a DW, which can induce spin 
scattering and affect current flowiii^iiLiSiiS. Third, the 
atomic structure (geometry) which affects both the elec- 
tronic and magnetic structures, and thus, the transmis- 
sion of these channels. Previous theoretical works present 
mutually conflicting results with methodologies that ei- 
ther used an oversimplified description of Ni electronic 
structureiSiSSiSi or idealized geometriesiSiSS. Here we 
present transport calculations across Ni nanocontacts de- 
scribing the electronic, magnetic, and atomic structure 
with ab-initio calculations^SiS^iS. Our results lead us to 
conclude that intrinsic ballistic MR is certainly not large 
in Ni nanocontacts. 

Spin- dependent transport formalism.- Transport 
through atomic-size metallic contacts is currently un- 
derstood in terms of elastic transport of non-interacting 
quasiparticles through a one-body potential that de- 
scribes their interaction with the constriction. In this 
approach, the conductance G is proportional to the 
quantum mechanical transmission T associated with 
the potential. On the other hand, the spontaneous 
breaking of the spin degeneracy in transition metal 
ferromagnets, which is due to electron-electron inter- 
actions, can also be properly understood in terms of 
a mean-field description, where quasiparticles interact 
with a spin-dependent self-consistent potential. Once 
the self-consistent field is determined for a given ge- 
ometry, the quantum-mechanical spin-dependent and 
energy-dependent transmission probabihty Xja'{E) can 
be obtained, and thereby the zero-bias conductance, 
using Landauer's formulaic.: 

G=^[T^^{EF) + T^l{E^) + T^liEF) + Ti^iEF)]. (1) 

In the above expression we only make explicit the de- 
pendence of %ja' on the spin channels, which we assume 
well-defined in the leads. 

Ab initio Cluster Embedded Calculations.- It is an ex- 
perimental fact that the chemical nature of the contact 
determines the conductance—. As a rule of thumb, the 
conductance of single-atom metallic nanocontacts can be 
as large as the number of valence orbitals, but, in prac- 



tice, is never larger than the number of valence elec- 
trons. A natural description of the problem is in terms 
of a localized atomic orbital basis, preferebly starting 
from first principles'^. In previous publications2i2^i2^ 
we have presented a method to pcrfom ab initio calcula- 
tions of quantum transport through atomic constrictions 
and molecules which is based on the code GAUSSIAN2&. 
Our approach has been successful in explaining experi- 
mental results in paramagnetic nanocontacts^iiSi. Here 
we take it a step further to study systems without spin 
degeneracy, like ferromagnetic nanocontacts. We solve 
the problem dividing the system in three different parts: 
left (L) and right (R) electrodes on one side and the con- 
tact region on the other. The spin-dependent one-body 
Hamiltonian is assumed fixed and homogeneous in the 
bulk electrodes, but it is determined self-consistently in 
the contact region subject to the appropriate magnetic 
boundary conditions. 

The density functional theory (DFT) calculations for 
the contact region are done with both LSDA and the hy- 
brid functional B3LYP'''. The LSDA results are robust 
against different basis sets and so we rely here on a mini- 
mal basis set with a core pseudopotential as described in 
previous work'^-'^. On the other hand, the B3LYP func- 
tional is more sensitive to the basis set due to its non- 
local exchange contribution. Therefore we employ here 
an all-electron basis set^. The electrodes are described 
by means of a semi-empirical tight-binding Bethe lattice 
model. With the appropriate parameters, the Bethe lat- 
tice can provide a geometry-independent description of 
the contacts with a bulk density of states (DOS) which 
is smoother than the real DOS and mimics an average 
over both disorder realizations and the actual electrode 
crystal orientations. Spin- mixing solutions are not con- 
sidered, i.e., Sz is a good quantum number. Thus, the 
last two terms in Eq. ^ do not give any contribution to 
the conductance. DW-like configurations are obtained 
for the adequate magnetic boundary conditions and the 
constraint Sz — 0- 

Results.- We restrict ourselves to the study of the last 
(first) plateau of conductance upon stretching (electrode- 
position) as, e.g., in the experiment of Viret et al.-^ (Sul- 
livan et alii). A reference atomic structure of the contact 
region has been initially taken like that shown in the in- 
set of Fig. n Following Viret et al^, we consider the 
narrowest (and most important) region to consist of two 
pyramids facing each other, formed along the (001) direc- 
tion, and with the two tip Ni atoms 2.6 A apart forming a 
dimer. Bulk atomic distances and perfect crystalline or- 
der are assumed otherwise. Ab initio simulations of the 
breaking process as the one shown in Fig. O support this 
choice. We stress that the section of the nanocontacts 
varies in the direction of the current flow. This is the 
situation in real nanocontacts and differs from perfect 1- 
dimensional systems, studied in Refs. |23, and from bulk 
systems studied by Van Hoof et ali^. In this regard, the 
geometries proposed by Bagrets et alm^ are closer to real 
nanocontacts, but are not backed up by experiments or 
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FIG. 1: (a) LSDA conductance per spin channel in the P 
configuration for the model nanocontact shown in the inset, 
(b) Same as in (a), but for the AP configuration. 



simulations. 

From the LSDA DOS projected on the tip atoms 
(not shown) we see that the sp orbitals are spin split- 
ted by less than 1 eV and that the minority (m) elec- 
trons are hybridized with the d levels, which are present 
at the Fermi energy. As a consequence, the LSDA 
DOS for the majority (M) electrons at the Fermi en- 
ergy is significantly smaller than that for the m ones. 
These results are compatible with LSDA first-principles 
calculations for systems with translational invariance— , 
have been properly accounted for in Anderson-like model 
Hamiltonians2£, but are in marked contrast to the Jsd 
model, usually invoked to understand large values of the 
MR in nanocontact s ^"^1 ^"^ . In these models the transmis- 
sion of itinerant s electrons is perfect in the ferromag- 
netic case while the d electrons do not contribute to the 
current since they are localized. The MR depends thus 
dramatically on the ratio between the spin splitting of 
the conducting s electrons and the Fermi energy when a 
DW is present. A large spin splitting needed to give large 
MR is, however, at odds with the actual Ni band struc- 
ture and this model must be ruled out from the outset 
to account for large MR in Ni nanocontacts. 

Figure (^ shows the LSDA conductance as a func- 
tion of energy for both up and down spin channels in 
two situations: (a) Parallel (P) and (b) antiparallel (AP) 
bulk magnetic arrangements. In both cases the self- 
consistent solution has been forced to respect the high 
symmetry of the nanocontact. In the AP case the self- 
consistent magnetization reverses abruptly between tip 
atoms. The resulting magnetic moment for the contact 
atoms is w 1.0/iB in both situations. This value is sig- 
nificantly larger than that obtained for bulk or surface 
atoms (« 0.6/^B and refiects the low coordination of the 
tip atoms forming the contact. In the P case the M 
channel is, for the most part, composed of a single sp 
orbital channel and conducts perfectly around the Fermi 
energy (set to zero) while the m channel is composed 
of three orbital channels (one sp- and two d-like, which 
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FIG. 2: (a) B3LYP conductance per spin channel in the P 
configuration for the model nanocontact shown in the inset of 
Fig. (b) Same as in (a), but for the AP configuration. 



conduct roughly the same), and exhibits a transmission 
strongly dependent on the scattering energy. In the AP 
case the system is invariant under the combined trans- 
formations that exchange L with R and | with J,, re- 
sulting in identical values for the conductance of the two 
spin channels, which now are composed of a dominant sp 
channel and a strongly diminished contribution of the d 
channels. The conductance ratio for this particular case 
is a; = 2.8/3.65 = 0.77. This yields MRi = 23% and 
MR2 = 30%, which is clearly below large MR claims^. 

LSDA provides a commonly accepted description of the 
electronic structure of bulk and surface ferromagnetism 
in transition metale^. However, the low coordination of 
atoms in nanocontacts might give rise to a further lo- 
calization of the d electrons (compared with bulk) and 
an increase of the magnetic moment. Since LSDA fails 
to describe properly localized electrons due to the self- 
interaction problem, various alternatives have been pro- 
posed to overcome this problem^, the most popular be- 
ing LDA-t-U2i, and SlOSi. We should point out that, 
while the results for the conductance reflect the DOS 
and look plausible, the fact that reported conductance 
histograms never show the lowest peak around 4e^//i at 
high magnetic fields'^'^ make us suspect that cither the 
chosen model for the atomic structure is not realistic or 
that the electronic structure given by LSDA, due to the 
problems mentioned above, does not provide the best ap- 
proximation for nanocontacts. We explore below both 
possibilities separately. 

An alternative approach to the electronic structure 
comes from the use of a hybrid functional like B3LYP 
which is a combination of Hartree-Fock and LSDA. The 
former is free from the self-interaction problem, but fails 
to include correlation, which is provided by the lat- 
ter. B3LYP happens to give a very good description of 
the electronic structure and local magnetic moments in 
NiO^ and LaaCuO^S^. With B3LYP the results for the 
conductance (see Fig. ^ are remarkably different in re- 
gard to the m channel. Now the d channels give a much 
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FIG. 3; Conductance of both spin channels for the P and 
AP configuration as a function of the stretching. The lines 
connecting the points are just to guide the eye. The insets 
show the relaxed contact geometry at different values of the 
displacement (OA, 1.5A, and SA from left to right). 



smaller contribution to the transmission at the Fermi en- 
ergy. In this case the MR is negative and its absolute 
value is even smaller (MRi = —11%) than the one ob- 
tained with LSDA. With B3LYP the bulk and surface 
magnetic moments are slightly higher than the LSDA 
ones while the magnetic moment for the tip atoms is 
roughly the same (« l/ie)- 

Since the m conductance evaluated at the LSDA level 
exhibits a strong dependence on the scattering energy, 
we study now whether or not different geometries can 
change the above results qualitatively. In an attempt 
to explore other realizations of the self-consistent poten- 
tial compatible with the magnetic boundary conditions 
and the experimental information, we perform ab initio 
structural relaxations as a function of the displacement 
between outer planes in the core cluster. To do so, we 
consider a cluster like that shown in Fig. Q] The inner 
atoms in the cluster (10 in total) are allowed to relax 
to local minimum energy configurations as we stretch. 
This results, logically, in lower energy solutions and in 
the loss of symmetry, so that the transmission in the AP 
case now becomes slightly spin-dependent. In Fig. |31the 
conductance at the Fermi energy per spin channel for the 
P and the AP configurations are shown as a function of 
the stretching up to the break-up point, starting from a 
slightly compressed nanocontact. From this figure we see 
that the conductance of the m channel for the P config- 
uration changes significantly upon small changes. The 
MR, on the contrary, barely changes as the nanocontact 
is stretched and is small, reaching vanishing values for 
the last points in Fig. |31 The conductance approaches a 
stable value around 2e^/h for both P and AP configura- 
tions. 

Discussion and conclusions.- As mentioned earlier, the 
maximum number of conducting channels in atomic-size 
contacts is roughly determined by the number of valence 
electrons of the contact atom(s). However, as shown 
above, this hypothetical upper limit is never reached. 



particularly for the m electrons, remaining essentially 
only one M channel and one m channel transmitting in 
the P case for stretched contacts. This result is impos- 
sible to predict without a full atomistic self-consistent 
calculation. The M channel is sp-type. Thus, this chan- 
nel transmits almost perfectly and evolves smoothly with 
the stretching of the contact giving a stable contribution 
7|f « 1 (see Fig. O The sp orbitals in the m channel 
are strongly hybridized with d-orbitals and, therefore, 
are more sensitive to the contact geometry. The con- 
tribution to the conductance of the latter, which form 
narrower bands, disappears with the stretching and dis- 
order, as expected (see Fig. I^J. On the other hand, in 
the AP configuration mostly one sp orbital channel per 
spin contributes. The conductance per spin channel lies 
thus in the vicinity of e^ /h, giving w 2e^ /h in total and 
is fairly stable during the last stage of the breaking of 
the nanocontact. 

To conclude, the reason behind the very small MR val- 
ues is the orbital (or geometric) blocking of most of the a 
priori available m channels in the P configuration due to 
the non-ideal geometry of the nanocontacts. The num- 
ber of bands at the Fermi energy in the case of a perfect 
mono-strand infinite ferromagnetic chainS^ is much larger 
than the number of non-zero eigenvalues of the transmis- 
sion matrix in a nanocontact (for a given basis and func- 
tional) . This phenomenon affects mainly the d bands (no 
DW involved) and, therefore, we call it orbital blocking. 
On the other hand, sp bands are less sensitive to geome- 
try. In fact, there can only be a significant and positive 
MR in the cases where the m channel in the P configura- 
tion conducts appreciably. This is the case for the highly 
symmetric nanocontact presented in Fig. ^within LSDA 
and, e.g., for the Ni chains studied in Refs. 22. In the P 
case, the number of M and m bands at the Fermi energy is 
1 and 6, respectively. When a DW is formed in the chain. 



5 out of 6 m channels are blocked, remaining 1 per spin, 
almost fully transmitting. As a result, one expects a ratio 
X ~ 2/7 resulting in MR2 ~ 250% for the ideal chainS^. 
The blocking of a number of channels is due to the fact 
that the d electrons, and not the sp electrons, are spin 
splitted. This could explain some early results^, but not 
recent ones^. Furthermore, to date, no evidence of chain 
formation in Ni has been reported. Even so, scatter- 
ing at the electrode-chain contact will always be present. 
For completeness, we have also performed calculations 
for Ni chains using the B3LYP functional. The number 
of bands crossing the Fermi energy is now reduced to 1-1-4 
compared to the LSDA results. This agrees with recent 
LSDA+U calculations reported by Wierzbowska et ali^tt 
where the two degenerate flat minority bands d^y, dx^-y^ 
are shifted downwards in energy because of the exchange 
interaction cancelling part of the self-interaction of the 
strongly localized electrons in these flat bands. However, 
this is not the reason for the drop in the m conductance 
seen in Fig|21 since the corresponding channel does not 
contribute to the conductance in the LSDA case either. 
In addition, non-coUincar DW's, not considered here, also 
reduce the MR^^. It is our believe that when observed, 
large values of the MR in Ni nanocontact might be due to 
magnetostriction effects and the corresponding formation 
of wider section contacts or to the presence of adsorbates, 
which modify the local electronic structure'^^ . 
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